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Abstract 

Pairing of oxygen holes into heavy bipolarons in the paramagnetic phase and their magnetic 
pair-breaking in the ferromagnetic phase [the so-called current-carrier density collapse (CCDC)] 
has accounted for the first-order ferromagnetic phase transition, colossal magnetoresistance (CMR), 
isotope effect, and pseudogap in doped manganites. Here we propose an explanation of the phase 
coexistence and describe the magnetization and resistivity of manganites near the ferromagnetic 
transition in the framework of CCDC. The present quantitative description of resistivity is obtained 
without any fitting parameters by using the experimental resistivities far away from the transition 
and the experimental magnetization, and essentially model independent. 
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Ferromagnetic oxides, in particular manganese perovskites, show very large magnetoresis- 
tance near the ferromagnetic transition. The effect observed in these materials was termed 
"colossal" magnetoresistance (CMR) to distinguish it from the giant magnetoresistance in 
metallic magnetic multilayers. The discovery raised expectations of a new generation of 
magnetic devices, and is a focus of extensive research aimed at describing the effect. Signifi- 
cant progress has been made in understanding the properties of CMR manganites, but many 
questions remain. The ferromagnetic metal-insulator transition in manganites has long been 
thought to be a consequence of the so-called double exchange mechanism (DEX), which re- 
sults in a varying bandwidth of electrons in the Mn 3+ d-shell as a function of temperature Q . 
More recently, it has been noticed P] that the effective spin interaction cannot alone account 
for CMR within the double-exchange model. In fact, there is strong experimental evidence 
for exceptionally strong electron-phonon interactions in doped manganites from the giant 
isotope effect 3[, the Arrhenius behaviour of the drift and Hall mobilities j^J in the para- 
magnetic phase above the Curie temperature, Tc, and other experiments. In view of this, 
Ref. jsj and some subsequent theoretical studies have combined DEX with the Jahn- Teller 
e-ph interaction with d-states arriving at the conclusion that the low-temperature ferromag- 
netic phase is a spin-polarised metal, while the paramagnetic high-temperature phase is a 
polaronic insulator. 

However, some low-temperature optical 0], electron-energy- loss (EELS) J(|, photoemis- 
sion [7] and thermoelectric js| measurements showed that the ferromagnetic phase of man- 
ganites is not a conventional metal. In particular, broad incoherent spectral features and 
a pseudogap in the excitation spectrum were observed. EELS confirmed that manganites 
were in fact charge-transfer doped insulators having p-holes as current carriers, rather than 
d(Mn 3+ ) electrons. Photoemission and x-ray absorption spectroscopies of Lai_ x Sr x Mn03 
also showed that the itinerant holes doped into LaMnOa have oxygen p-character. Further, 
CMR has been observed in the ferromagnetic pyrochlore manganite T^M^Oy j^, which 
has neither the mixed valence for DEX magnetic interaction nor the Jahn- Teller cations 
such as Mn 3+ . 

These and other observations Q], in particular the fact that some samples of ferromag- 
netic manganites manifest an insulating-like optical conductivity at all temperatures 11], 
clearly rule out DEX as the mechanism of CMR. The earlier of the above observations 
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led to a novel theory of ferromagnetic/paramagnetic phase transition and 
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FIG. 1: Bipolaron model of CMR: pairs (BP) are localised on impurity levels in the paramagnetic 
phase, where the only current carriers are single thermally excited polarons. If the exchange 
interaction JSa between p-hole polarons and ordered manganese spins exceeds the pair binding 
energy A , the pairs break at T < Tq because the spin-up polaron sub-band sinks abruptly below 
the bipolaron level. The ferromagnetic state is a polaronic conductor. 



CMR, based on the so-called current-carrier density collapse (CCDC) |l2j, confirmed by 
later observations. In CCDC model the p-holes are bound into heavy bipolarons above the 
Curie temperature Tq due to the Frohlich electron-phonon interaction, which is written in 
the real-space representation as 



nnV=T,| 



(1) 



where £ n ' is the ion displacement operator, and the form of electron-phonon interaction is 
specified via the force function |l3j / n /(n). The latter is defined as the force with which an 
electron in state |n) interacts with the ion degree of freedom £ n /. 

The resistivity peak and CMR are the result of the magnetic pair-breaking below Tc, 
Fig. 1, caused by the p — d spin-exchange interaction, J p d, described as 

H pd = -(2N)- 1 J P dS z m (ci^ - cl lCnl ). (2) 

n,m 

Here is the z-component of Mn 3+ spin on site m, c n | and c n i annihilate a p(|, |)-hole 



on the oxygen site, n, with spin up and down, respectively, and iV is the total number of 
unite cells. 

Different from cuprates, hole on-site or more extended intersite oxygen bipolarons are 
much heavier in manganites because the e-ph Frohlich interaction, Eq.(l), is stronger ^| 
and the band structure is less anisotropic. They are readily localised by disorder, so it is 
mainly thermally-excited single polarons that conduct in the paramagnetic phase. Upon 
temperature lowering single polarons polarize manganese spins at Tc via the exchange in- 
teraction J p d, and the spin polarization of manganese ions breaks the bipolaronic singlets 
creating a spin-polarized polaronic conductor. 

The CCDC model has explained CMR in the experimental range of external magnetic 
fields HQ. More recently, the theory has been further confirmed experimentally. In 
particular, the oxygen isotope effect has been observed in the low-temperature resistivity of 
La .75Ca .25MnO3 and Ndn 7 Sr .3MnO3 and explained by CCDC with polaronic carriers in 
the ferromagnetic phase Q. The current-carrier density cohapse has been directly observed 
using the Hall data in Lao.67Cao.33Mn03 and Lao.67Sr .33Mn03 |17l |. and the first order phase 
transition at Tc, predicted by the theory 3], has been firmly established in the specific 
heat measurements Q|. Importantly the character of the magnetic phase transition in 
Tl2iVIn207 pyrochlores has also been determined to be the first order 19( and attributed to 
the tendency of small polarons to phase separation at finite carrier density. Indeed recent 
Monte Carlo simulations j^j of lattice polarons with anisotropic e-ph interactions and the 
realistic long-range Coulomb repulsion show diverse mesoscopic textures in the adiabatic 
limit, where spatially disordered pairs (i.e. bipolarons) dominate at finite doping. 

On the other hand, resistivity and the magnetization ot some La ,Ca . 3 Mn 1 _ J Ti J O 3 sam- 
pies showed a more gradual (second-order like) transition |21|. Also the coexistence of ferro- 
magnetic and paramagnetic phases near the Curie temperature observed in tunnelling 
and other experiments has not yet been addressed in the framework of CCDC. Here we 
show that the diagonal disorder, which is inevitable with doping in those solid-state solu- 
tions, explains both the phase coexistence and the resistivity/magnetization behaviour near 
the transition. 

The mean-field equations 



12j describing p-hole polaron atomic density, n, polaron, m, 
and manganese, a, reduced magnetizations, and the chemical potential fi = ksTlny are 
readily generalized taking into account a random distribution of the bipolaron binding energy 
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FIG. 2: Experimental magnetizations in Lao.7Cao.3Mii2_ ( $Ti ( $03 (symbols 21 j) compared with 
Eq.(4) (lines) with T = 8K, T c = 278K for 5 = 0.00; T = 8K, T c = 264K for S = 0.01; and 
T = 18K, T c = 224K for 5 = 0.03. 



5 = A/(2J pd ) across the sample, 



m = 6ycosh((Ji/t), 
rrii = rii tanh((7j/t), 
Oi = B 2 (m,i/2t), 

2/ 2 = ^exp(-2^A), (3) 



where t = ksT '/ J p d is the reduced temperature, B$ is the Brillouin function, x is the number 
of holes at zero temperature in p-orbital states, which are 3-fold degenerate. The subscript 



121 with different Curie 



i means different parts of the sample with different Si and hence 
temperatures, Tci, owing to disorder. 

While averaging these simple equations over a random distribution of 5i is rather cum- 
bersome, one can apply a simplified approach using the fact that the phase transition in a 
homogeneous system is of the first order in a wide range of 8 [l^] . Taking o"j ~ O (Tci — T) 
and rii ps xQ(T C i — T) + v / 2xexp(— A/(2k B T))Q(T — T C i), and averaging both quantities 
with the Gaussian distribution of random TciS around the experimental Tc we obtain an 
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averaged manganese magnetization 

a(T) = ^erfc ) , (4) 

where A is the average bipolaron binding energy, Q(y) = 1 for y > and zero for y < 
, and erfc(;z) = J°° dy exp(— y 2 ) . The CCDC with disorder, Eq. (4) fits nicely the 



experimental magnetizations 2l( near the transition with physically reasonable T of the 
order of 10K, depending on doping, Fig. 2. Hence, we believe that the random distribution 
of transition temperatures with the width T across the sample caused by the randomness of 
the bipolaron binding energy is responsible for the phase coexistence near the transition as 
seen in the tunnelling experiments (3]. We note that some dro p o f magnetization at low 
temperatures as seen in Fig. 2 might be caused by domain walls 

Resistivity of inhomogeneous two-phase systems has to be calculated numerically. Nev- 
ertheless, the comprehensive numerical simulations are consistent with a simple analytical 
expression for the resistivity of the binary mixture, 

P = fffi, (5) 



which is valid in a wide range of the ratios pi/p2 [24]]. Here pi.2 is the resistivity of each phase, 
respectively, and v is the volume fraction of the second phase. The expression Eq.(5) is a 
homogeneous function of p\$, satisfying the duality relation. Scaling arguments [2J| proved 
that the expression is exact near a percolation threshold in two dimensions. Numerical 



analysis 24] has shown that Eq.(5) describes the effective resistivity of 2D random systems 
even far away from the percolation threshold, if the resistivity ratio p\j pi is not extremely 
large (< 20). The same expression is also asymptotically correct for any 3D system, if 
|1 — Pxj pi\ "C 1, and for specific 3D lattice structures [3], if p\j P2 ^> 1- Of course, there 
is no universal formula for any material with randomly distributed phases. Generally, one 
could write f(p) = (1 — v) f{p-\)+yf{p<>), where f(x) is a model function (for a comprehensive 



list of mixture formulae see Ref. 
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281]). Equation (5) corresponds to the average of lnp in 
isotropic mixtures providing a qualitatively reasonable and numerically accurate description 
of the effective resistivity in many physically important cases (see below). 

In the framework of CCDC, the resistivity of the paramagnetic phase is pi(T) = 
f(T) exp(A/2fceT) and the resistivity of the ferromagnetic phase is pi(T) = <f)(T), where 
f(T) and 4>(T) are polynomial functions of temperature depending on the scattering mech- 
anisms. Well below the transition <f>(T) can be parameterized as 4>{T) = p + aT 2 , and 
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FIG. 3: CCDC model (Eq.(5), lines) describes the experimental resistivity near the ferromagnetic 
transition in Lao.7Cao.3Mn2_<5Ti,503 (squares 2^|), if the phase coexistence caused by disorder is 
taken into account. No fitting parameters are used in Eq.(5) but the experimental resistivity well 
below and well above the transition and the experimental magnetization. 



f(T) = bT well above the transition, where the temperature independent parameters po, a, 
A/2 and b are taken directly from the experiment The microscopic ori gin of po, a, b 
and alternative parametrization formulae have been discussed e.g. in Ref. QH), and are 
not an issue here. The volume fraction v of the ferromagnetic phase is simp 



y the relative 



magnetization in our model, v = cr(T), also available from the experiment 2l|. As a result, 
Eq. (5) provides the quantitative description of p(T) in the transition region without any 
fitting parameters by using the experimental resistivities far away from the transition and 
the experimental magnetization j2l|, as shown in Fig. 3. On the other hand, fitting the 
ferromagnetic-phase resistivity with a magnetic scattering (p™ 5 T 4,5 ) leads to an unrealistic 
doping- dependent coefficient p™ 5 that is changing with doping by more than five orders of 
magnitude (see the Table in Ref. |2 1] ) . Note, that if one were using an estimate, p oc 1/n, 
where n is the average single-polaron density, one would obtain 



p 1 (T) oc erfc 



T-T, 



c 



+ (2/x) 1 / 2 e- A / 2fcBT erfc 



T c -T 



(6) 



r j w ' v r 

This expression can also fit the experimental curves, but with a value of Tc, which turns out 



to be smaller than that in the magnetization, Eq. (4), by several tens of degrees Kelvin [30]. 
The latter expression corresponds to a linear expansion of Eq. (5) in powers of 1 — p\j pi- It 
is easy to see why Eq. (5), when compared with Eq. (6), resolves the problem of different 
Tcs in the magnetization and resistivity, thus providing a parameter-free description of 
experimental p(T). If we take p\ ^> p 2 , the resistivity at the magnetic transition (i.e. for 
v = a = 1/2), p = P2A/P1/P2, calculated with Eq. (5) turns out larger than the resistivity 
p ~ 2p 2 calculated with the perturbation expression, Eq. (6). It is important that the present 
description of resistivity does not depend on a particular model, but on the assumption that 
the ferromagnetic transition is the first order. The CCDC has provided a basis for this 
assumption in terms of the microscopic model. 

In summary, we have shown that the conventional DEX model, proposed half a century 
ago and generalized more recently to include the electron-phonon interaction, is in conflict 
with a number of recent experiments. Among these experiments are the site-selective spec- 
troscopies, which have shown unambiguously that oxygen p-holes are the current carriers 
rather than d-electrons in ferromagnetic manganites. Also, some samples of ferromagnetic 
manganites manifest an insulating-like optical conductivity at all temperatures, contradict- 
ing the DEX notion that their ferromagnetic phase is metallic. On the other hand, the 
pairing of oxygen holes into heavy bipolarons in the paramagnetic phase and their magnetic 
break-up in the ferromagnetic phase has explained the colossal magnetoresistance, the iso- 
tope effects, and the pseudogaps observed in doped manganites. It also explains the CMR 
in systems where DEX simply cannot exist, like manganese pyrochlores The CCDC 
theory of CMR predicts the first-order phase transition and allows the present simple ex- 
planation of the coexistence of high and low-resistive phases. It explains the temperature 
dependencies of the magnetization and the resistivity near the transition as the result of the 
unavoidable disorder and transport through the two-phase mixture in doped manganites. 
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